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We present a quantitative shadowgraphic method which can measure the density of a laser-generated plasma 
in air with sensitivity and resolution comparable or better than traditional interferometric techniques. Simul- 
taneous comparison of both shadowgraphy and interferometry has been carried out allowing the experimental 
evaluation of the reliability of the shadowgraphic method. © 2008 Optical Society of America 
OCIS codes: 100.5070, 350.5400. 
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Many branches of physics face the problem of the accu- 
rate characterization of small distributions (d ~ 20fim) 
of low density plasmas (n e < 10 19 cm -3 ), such as laser- 
produced plasma channels in dielectrics, 1 and coronal 
plasma jets generated in pulsed power devices (e.g. X- 
pinch 2 ) or during high-intensity laser-solid interactions. 3 
Optical interferometry is the most favoured method 
for characterizing under-critical plasmas. 4, 5 However, to 
probe a plasma slab of 20/im of thickness and density 
n e = 10 19 cm -3 with visible light, a set-up with sen- 
sitivity far better than ~ A/20 is necessary, which is 
usually quite a demanding requirement. Moreover, many 
phase retrieval methods can not exploit the full resolu- 
tion of the imaging system, the smallest resolvable de- 
tail is limited by the size of the interference fringes. Re- 
cently, the plasma produced in laser-induced filaments 
in transparent solids has been characterized by means of 
measurements of the inverse bremsstrahlung absorption 6 
of a probe pulse propagating perpendicularly to the fil- 
ament. 7,8 However, the density measurements obtained 
by these techniques heavily rely either on assumptions 
on the electron-ion collision time, 7 or on an independent 
interferometric measurement, 8 which is again limited in 
sensitivity and spatial resolution. 

In this Letter, we present a set-up based on quanti- 
tative refractive shadowgraphy 9-11 which can overcome 
these limitations and achieve a remarkable spatial reso- 
lution (~ 25/im) and sensitivity (An e ~ 5 • 10 17 cm -3 ), 
useful for the characterization of laser-induced plasmas 
in air. 

According to geometrical optics in a homogeneous 
medium, optical rays belonging to a collimated unifrom 
beam of intensity Iq propagate parallel to each other in 
the longitudinal direction z. The rays hitting a transpar- 
ent, optically thin object are deflected by its transverse 
gradients of the refrcctivc index and result in a trans- 
verse intensity pattern I(x, y) (the shadowgram)at an 
obsevation plane placed at a distance L from the object. 
By assuming negligible diffraction effects, small angle de- 
flection, and the thickness of the sample Az <C L, it is 
possible to relate the variation of the transverse optical 



thickness of the sample to the shadowgram 
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where n(x, y, z) is the refractive index distribution of the 
phase object. It is easy to recognize in Eq. 1 the Poisson 
equation: the unknown is the optical thickness of the ob- 
ject, while the source term is proportional to the normal- 
ized intensity variation of the shadowgram respect to the 
intensity of the incoming plane wave. For radially sym- 
metric samples, the optical thickness distribution can be 
inverted by means of the Abel inversion transformation 12 
resulting in the actual radial distribution of refractive 
index n(r), which in turn is proportional to the plasma 
density radial profile. 

Theoretically, the accuracy of the shadowgraphic 
method is limited only by the impact of diffraction in 
free space propagation. 13 This can be minimized by tak- 
ing L small enough or by using a short-wavelength light 
source as object backlighter. From the practical view- 
point, two main problems have to be solved in order to 
achieve a satisfactory accuracy, namely: i) reduction of 
large-scale correlated noise sources (e.g. air turbulence), 
and ii) coding of fast numerical schemes for the resolu- 
tion of the Poisson equation on large arrays of data. 

To test the reliability of the shadowgraphic method, 
we used the set-up shown in Fig. 1. A pulse from our 
amplifed sub-ns laser source (EKSPLA SL312, deliver- 
ing 150 ps pulses at a wavelength of 1064 nm and energies 
up to 250 mJ) is expanded and divided into two channels 
by a semi-reflecting dielectric mirror. The reflected beam 
(pump) is focused in air using a / = 300 mm lens to ex- 
cite a line of plasma bubbles through laser breakdown. 
The transmitted beam (probe) is frequency doubled us- 
ing a 1 cm thick KDP crystal and passed through an 
adjustable delay line (delay from to 3 ns respect to 
the time of plasma bubble formation) and a spatial filter 
before probing orthogonally the plasma bubbles. The fre- 
quency conversion of the probe was required in order to 
reduce the impact of the inverse bresmstrahlung absorp- 
tion. The probe is then collected by a / = +50 mm lens 
and divided in two channels by a 50/50 beamsplitter. 
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On one channel, the shadowgram of a plane at distance 
L = 1.5 mm from the bubbles is recorded. On the sec- 
ond channel a Nomarsky interferometer was built using 
0.5 degree split angle Wollaston prism. The camera is 
focused on a plane containing the line of plasma bubbles 
to minimize distorsions. 4 Both shadowgram and interfer- 
ogram were recorded simultaneously using a dual frame 
grabbing software and by mutually synchronizing the 
cameras. Each measurement consisted of taking a frame 
of the plasma bubbles followed by a measurement of the 
background taken while stopping the pump beam. The 
shadowgraphic images were processed using a specifically 
developed fast Poisson solver using optimized successive 
over relaxation method, while the interferograms were 
processed by means of the IDEA software. 14 The whole 
set-up was enclosed in a box to avoid any disturbance in 
the air flow which could affect the shadowgraphic data. 
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Fig. 2. a) Shadowgram of a chain of expanding plasma 
bubbles excited by a focused laser pulse in air (to = 150 
ps, A = 1064 nm, E = 10 mJ). Time delay of the probe 
pulse (to = 110 ps, A = 532 nm): +2.0 ns. Image taken at 
L = 1.5 mm from the object plane, b) and c): Retrieved 
phase maps resulting from the shadowgraphic method 
and from the interferometric set-up, respectively. The 
double image in the interferometric map is a peculiarity 
of the modified Nomarski interferometer we used. The 
numbers refer to the bubbles displayed in Fig. 3. 



Fig. 1. Sketch of the experimental set-up 

A typical outcome of the processing of the data ac- 
quired by the two cameras is illustrated in Fig. 2. A 
qualitative comparison between the interferometric and 
shadowgraphic measurements (Fig. 2.b and 2.c) reveals 
that the shadowgraphic method allows an accurate map- 
ping of the array of laser-generated plasma bubbles in 
greater detail compared to the interferometric method. 

Cross-section optical thickness profiles show that 
the shadowgraphic retrieval matches the interferometric 
measurement also on a quantitative basis (phase profiles, 
Fig. 3). For low to medium density plasma bubbles (Fig. 
3. a and 3.b) the accuracy of the shadwographic retrieval 
is remarkably good. The peak values measured with both 
methods are the same within the measured large-scale 
fluctuation of the baselines (±0.1 rad). Regarding the 
spatial width of the peaks, they match within the esti- 
mated resolution of the interferometric set-up (1 fringe 
separation = 28/im). Abel inverted plasma density pro- 
files (Fig. 3.d and 3.e) tend to enhance the differences 
near the peak value and reduce the impact of irregu- 
lar baselines, however the agreement between the two 
methods is very good. The discrepancy found for low- 
density profiles has to be attributed mainly to the lower 
resolution of the interferometric set-up. For higher den- 



sity profiles (Fig. 3.c), the shadowgraphic method appar- 
ently fails in reproducing accurately the interferometric 
phase data. In fact, the outcome of shadowgraphy shows 
the plasma peak sitting on a large-scale parabolic back- 
ground affecting the absolute peak value of the phase 
shift. Notice however that the peak is indeed reproduced 
in great detail, as confirmed by the Abel transformed 
data (Fig. 3.f). 

The appereance of a parabolic backgound in the shad- 
owgraphic phase data can be explained quantitatively by 
assuming significant absorption of the probe beam in the 
plasma bubble due to inverse bremssthrahmg effect. 6 In 
fact, by suitable Fourier filtering of the interferograms 
we were able to reconstruct a low resolution L = shad- 
owgram which clearly shows an absorption maximum of 
12 ± 5% in the position of peak no. 3. We simulated the 
impact of absorption by numerically generating shadow- 
grams of a Gaussian distribution of plasma in the frame 
of a paraxial beam propagation model. If the absorption 
losses are not taken into account in the background frame 
(i.e. Io(x,y) =const.) the retrieved phase map will de- 
velop a parabolic background such as the one observed 
experimentally. The best fit of the experimental data 
(not shown) is obtained for a peak absoprtion of 24%. If 
the data are processed by using as a reference a shadow- 
gram at L = 0, the parabolic background can be removed 
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Fig. 3. Quantitative comparison of the shadowgraphic 
and interferometric method. a)-c) Phase profiles along 
the vertical direction of bubbles f, 2 and 3 of Fig. 2, 
respectively. d)-f) Retrieved electron density profiles of 
the same bubbles. 



almost completely (see Fig. 3.c). 

The beam propagation model was also used to evalu- 
ate the ultimate sensitivity and resolution of the shad- 
owgraphic method. The simulations 15 basically confirm 
previous results found in the literature 9, 13 which predict 
a negligible impact of diffraction as long as i) the ob- 
servation distance L is much smaller than the absolute 
value of the effective focal length of the phase object, and 
ii) the far-field diffraction spot is much smaller than the 
transverse dimensions of the object. These requirements 
are in contrast to the requirement of large L to achieve 
high sensitivity therefore, a tradeoff between resolution 
and sensitivity has to be found. For a green light probe 
recorded with an 8-bit CCD camera, our simulations in- 
dicate that L between 1.0 and 1.5 mm is a good com- 
promise for the accurate retrieval of plasma bubbles with 
density in the range 0.5-20 10 18 cm -3 and FWHM diam- 
eters greater than ~26-33 (im (accuracy better than 10% 
in both peak phase and width retrieval). The ultimate 
resolution of the method can be estimated by equating 
the spot size of the far-field diffraction to the actual size 
of the object. For L in the chosen range, this gives a 
value of FWHM of the gaussian peak in between 20 and 
26 /im. 

In conclusion, we have shown that quantitative refrac- 



tive shadowgraphy is a reliable and useful technique for 
the accurate measurement of low-density, undercritical, 
radially symmetric plasma distributions. Compared to 
the existing literature, the novelty of our method re- 
sides on the possibility to retreive general two dimen- 
sional phase maps without the constraint of periodicity 
of the pattern. 11 We finally point out that, in contrast 
to shear intcrfcromctry, the method could in principle 
be scaled to probe pulses of a few optical cycles, thus 
opening interesting perspectives for the measurement of 
time resolved plasma dynamics during femtosecond pulse 
filamentation in dielectric media. 1 
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